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Abstract. The cosmic energy density in the form of radiation before and during Big 
Bang Nucleosynthesis (BBN) is typically parameterized in terms of the effective number 
of neutrinos iV e gf. This quantity, in case of no extra degrees of freedom, depends upon the 
chemical potential and the temperature characterizing the three active neutrino distributions, 
as well as by their possible non-thermal features. In the present analysis we determine the 
upper bounds that BBN places on N e g from primordial neutrino-antineutrino asymmetries, 
with a careful treatment of the dynamics of neutrino oscillations. We consider quite a wide 
range for the total lepton number in the neutrino sector, r) v = r) Ve + r) v + t] Vt and the 
initial electron neutrino asymmetry r?™, solving the corresponding kinetic equations which 
rule the dynamics of neutrino (antineutrino) distributions in phase space due to collisions, 
pair processes and flavor oscillations. New bounds on both the total lepton number in the 
neutrino sector and the v e — v e asymmetry at the onset of BBN are obtained fully exploiting 
the time evolution of neutrino distributions, as well as the most recent determinations of 
primordial 2 H/H density ratio and 4 He mass fraction. Note that taking the baryon fraction 
as measured by WMAP, the 2 H/H abundance plays a relevant role in constraining the allowed 
regions in the r\ v — 77™ plane. These bounds fix the maximum contribution of neutrinos with 
primordial asymmetries to iV e fj as a function of the mixing parameter #13, and point out the 
upper bound N e g < 3.4. Comparing these results with the forthcoming measurement of N e g 
by the Planck satellite will likely provide insight on the nature of the radiation content of 
the universe. 
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1 Introduction 



The dynamics of neutrino oscillations in the early universe has been extensively studied in 
the scientific literature and is expected to have produced an efficient mixing of flavor neutrino 
distributions at a temperature T 7 ~ 1 MeV, thus shortly before the freezing of the neutron-to- 
proton density ratio and the onset of Big Bang Nucleosynthesis (BBN). Indeed, while flavor 
neutrino conversions are suppressed by matter effects at larger temperatures, at T 7 ~ 10 
MeV the atmospheric mass-squared difference Am^ tm triggers efficient Vn — v T mixing, as 
well as v x — v e (x = fi, r) conversions if the angle #13 is not vanishing and sufficiently large. 
Finally, at T 7 < 3 MeV oscillations driven by Am^ set on and are large enough to achieve 
strong flavor conversions before BBN [1-4]. 

These results have a major impact on the possible values for the cosmological lepton 
asymmetry stored in each neutrino flavor, which in analogy with the baryon-antibaryon 
asymmetry parameter r]b = (rib — n b)/ n ~ii can be parameterized by the number density ratios 



Based on the equilibration of lepton and baryon asymmetries by sphalerons in the very early 
universe, such a neutrino asymmetry should be of the same order of the cosmological baryon 
number rjb = 273.93 Q,bh? 1CP 10 , which is restricted to be a few times 1CP 10 by present obser- 
vations, such as 7-year data from the WMAP satellite and other cosmological measurements 
[5]. Nevertheless, a cosmological neutrino asymmetry orders of magnitude larger than this 
value is still an open possibility, with implications on fundamental physics in the early uni- 
verse, such as their potential relation with the cosmological magnetic fields at large scales [6]. 
In particular, a non-zero lepton asymmetry leads to an enhanced contribution of neutrinos 
to the energy density in the form of radiation p r , which, after the e + e~ annihilation phase, is 
usually parameterized as p r / Pl = 1 + 7/8(4/ll) 4 / 3 iV eff . The parameter N eS is the "effective 
number of neutrinos" whose standard value is 3 in the limit of instantaneous neutrino de- 
coupling. Interestingly, recent data on the anisotropies of the cosmic microwave background 
from WMAP [5] and the primordial 4 He abundance [7-9] seem to favor a value of N e g > 3, 
although with large errorbars. 

If neutrinos are in kinetic and chemical equilibrium, their distribution of momenta is 
parameterized by a temperature and a well defined chemical potential fj, Ua , and each flavor 
neutrino asymmetry can be expressed in terms of the corresponding degeneracy parameter 
iv a = liv a /T Ua as 



In this case, if flavor oscillations enforce the condition that all £„ a are almost the same 
during BBN, the stringent bound on the electron neutrino degeneracy, which directly enters 
the neutron/proton chemical equilibrium [10], applies to all flavors. The common value 
of the neutrino degeneracies is restricted to the range —0.021 < £ Va < 0.005, a = e, [A, r 
[11] (see also [12-19] for other analyses). This in turn also implies that, if neutrinos indeed, 
reach "perfect kinetic and chemical equilibrium before they decouple, any large excess in cosmic 
radiation density, if observed, must be ascribed to extra relativistic degrees of freedom since 
the additional contribution to radiation density due to non vanishing £„ a is very small. 
Assuming that the neutrino distributions are given by their equilibrium form, the BBN bound 
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on the neutrino degeneracies leads to the following upper limit on the excess contribution to 



which is tiny, even compared with the value of N e g = 3.046 found solving the neutrino kinetic 
equations in absence of asymmetries [20] . 

Neutrino kinetic and chemical equilibrium is maintained in the early universe by purely 
leptonic weak processes such as neutrino-neutrino interactions, v-e^ scatterings and pair 
processes, vv -H- e + e~, whose rates become of the order of the Hubble parameter at T 7 ~ 
2 — 3 MeV. Baryons play no role in this concern due to their very low density. At lower 
temperatures weak interactions are no longer effective, neutrinos decouple from the rest of 
the primeval plasma and their distribution in phase space is frozen out. Flavor oscillations 
driven by Am^ tm take place when neutrinos are still fastly scattering off the surrounding 
medium, so that the changes in their distribution due to oscillations are efficiently readjusted 
into an equilibrium Fermi-Dirac function. Instead, flavor conversions due to A?rig ol and 
#12 occur around neutrino decoupling. This implies, at least in principle, that if neutrinos 
succeed in achieving comparable asymmetries in all flavors before BBN, their distributions 
might acquire distortions with respect to equilibrium values due to inefficient interactions. 

This can be easily understood by a simple example. Suppose that at temperatures 
higher than 2 — 3 MeV we start with a vanishing total asymmetry, but = — 2r/Jf 7^ and 
we artificially switch-off scattering and pair processes. Due to solar-scale oscillations, in the 
case of maximal mixing asymmetries in each flavor will eventually vanish, but the neutrino 
distributions will not correspond to equilibrium, since averaging two equilibrium distributions 
with a different chemical potential does not correspond to a Fermi-Dirac function. Only 
scatterings and pair-processes can turn it into an equilibrium distribution with, in this case, 
zero chemical potential. 

Though this example is quite extreme and unrealistic, nevertheless, it tells us that the 
interplay of neutrino freeze-out and Am^ ol oscillation phases might deserve a more careful 
scrutiny, as first discussed in [21]. In fact, depending on the initial flavor neutrino asymme- 
tries and the value of #13, the final neutrino distributions at the onset of BBN might show 
non-thermal distortions which change the neutron-proton chemical equilibrium due to the 
direct role played by electron (anti)neutrinos. Moreover, this corresponds to an asymmetry- 
depending parameter N e s > 3 which is due to inefficient entropy transfer to the electro- 
magnetic plasma and is not given by the equilibrium value of eq. (1.3). As pointed out in 
[21] these features will prove to be quite important in case of large initial asymmetries and 
opposite values for v e and v„ T chemical potentials. It was shown that with fine-tuned initial 
asymmetries the BBN bound could be respected and at the same time an excess radiation 
density could survive, corresponding to values of AN e g of order unity or larger. 

In the present work we extend the analysis of [21] in two ways. First, we consider 
a wider range of values of the initial neutrino asymmetries and solve their evolution with 
the corresponding kinetic equations, including both collisions and oscillations. Moreover, 
the obtained shape of the neutrino distributions is then plugged into the BBN dynamics 
allowing, by the comparison between the theoretical results and the experimental data on 
primordial abundances of deuterium and 4 He, to find more accurate bounds on the total 
lepton asymmetry stored in the neutrino sector, as well as the way it was distributed at some 
early stage in the v e and v x flavors. In fact, for all initial values of neutrino asymmetries which 
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are compatible with BBN bounds, the distortions in the neutrino distribution are typically 
quite small, see Section 2, so that it is accurate enough for our purposes to parameterize them 
in terms of a Fermi-Dirac function with two time-dependent parameters, which correspond 
to the first two moments of the actual distribution: an effective chemical potential £„ a and 
an effective temperature T Uct , or equivalently, the asymmetry in each flavor and the energy 
density. 

The paper is organized as follows. After setting in Section 2 the formalism of kinetic 
equations which rule the evolution of neutrino distributions written in the standard density 
matrix formalism and showing an example of their dynamics, we then study the BBN con- 
straints on neutrino asymmetries in Section 3. In particular, we discuss the experimental 
data which are used in our analysis for 2 H/H and 4 He mass fraction Y p , as well as the way 
we have modified the public BBN numerical code PArthENoPE [22-24] to track neutrino evo- 
lution. Finally, we report the bounds on initial (at T 7 ~ 10 MeV) neutrino asymmetries or 
their final values after flavor oscillation phase. In Section 4 we give our concluding remarks. 

2 The dynamics of neutrinos in the early universe with primordial asym- 
metries 

Our first aim is to calculate the evolution of the three active neutrino distributions in the 
epoch of the universe right before BBN, when these particles were interacting among them- 
selves and with electrons and positrons. The corresponding weak collision rate decreases very 
fast with the expansion until neutrinos decouple at T 7 ~ 1 MeV. At the same time, it was 
shown in [1-4, 25] that for the neutrino mixing parameters currently allowed, flavor oscilla- 
tions start to be effective at similar temperatures. In such a case the best way to describe 
neutrino distributions is to use matrices in flavor space [26, 27]. Since we consider only the 
three active neutrino species, we will need 3x3 matrices in flavor space g p for each neutrino 
momentum p, where the diagonal elements are the usual flavor distribution functions (oc- 
cupation numbers) and the off-diagonal ones encode phase information and vanish for zero 
mixing. 

Oscillations in flavor space of the three active neutrinos are driven by two mass-squared 
differences and three mixing angles bounded by the experimental observations in the following 
ranges: the "solar" Am s 2 ol = 7.59±Jj;$ x 10" 5 eV 2 , the "atmospheric" |Am 2 tm | = 2.40+^2 x 
1CP 3 eV 2 and correspondingly, the large mixing angles sin 2 0\2 = 0.32^q'q3 and sin 2 $23 = 
0.50^g'} 3 (2a ranges from [28]). On the other hand, the third angle is quite small, sin 2 #13 < 
0.053 (3a) [28], even compatible with a vanishing value, though a mild evidence for sin 2 #13 > 
has been found in global data analyses (see e.g. [29, 30]). 

The equations of motion (EOMs) for g p are the same as those considered in reference 

[21], 

i^ = [n p ,Q P ]+C[Q p ,Q p ], (2.1) 

and similar for the antineutrino matrices g p . The first term on the r.h.s. describes flavor 
oscillations, 

n p = ^ + V2G F (-^E + e -,), (2.2) 

where p = |p| and M is the neutrino mass matrix (opposite sign for antineutrinos) . Matter 
effects are included via the term proportional to the Fermi constant G-p, where E is the 3x3 
flavor matrix of charged-lepton energy densities [26] . For our range of temperatures we only 
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need to include the contribution of electrons and positrons. Finally, the last term arises 
from neutrino-neutrino interactions and is proportional to g — g, where g = j £ p d 3 p/(27r) 3 
and similar for antineutrinos. For the relevant values of neutrino asymmetries this matter 
term dominates and leads to synchronized oscillations [2-4]. The last term in eq. (2.1) 
corresponds to the effect of neutrino collisions, i.e. interactions with exchange of momenta. 
Here we follow the same considerations of ref. [21], where the reader can find more details 
on the approximations made and related references. In short, the collision terms for the off- 
diagonal components of g p in the weak-interaction basis are momentum-dependent damping 
factors, while collisions and pair processes for the diagonal g p elements are implemented 
without approximations solving numerically the collision integrals as in [20]. These last 
terms are crucial for modifying the neutrino distributions to achieve equilibrium with e 
and, indirectly, with photons. 

We have solved numerically the EOMs for the matrices in flavor space of neutrinos 
and antineutrinos with non-zero initial asymmetries. The expansion of the universe is taken 
into account using comoving variables as in [2], where it was shown that flavor oscillations 
between muon and tau neutrinos take place at T 7 > 10 MeV, when interactions are very 
effective. For any initial values of the muon or tau neutrino asymmetries, the combined 
effect of oscillations and collisions is able to equilibrate the two flavors and hence leads to 
r/Jf = 77™ = 77™ . Therefore, we start our numerical calculations for each case at T = 10 MeV 
and initial parameters 77™ and the total asymmetry j\ v = 77™ + 2 r/Jf . Note that hereafter by 
r\ v we denote the initial value of the total asymmetry, which is kept constant until the onset 
of e + — e~ annihilations, when it is diluted by the increase of the photon density (as in the 
case of the baryon asymmetry 77;,). 

All neutrino mixing parameters, except for #13, are taken as the best-fit values in [28]. 
Modifying these parameters within their allowed regions does not affect our results. Instead, 
the value of #13 plays an important role in the evolution of the neutrino asymmetries [21]. 
We thus consider either #13 = or sin 2 #13 = 0.04, a value close to the upper bound from 
neutrino experiments [28-30]. 

Let us describe the evolution of the flavor asymmetries with a specific example. As 
shown in [21], large initial values of flavor primordial asymmetries could satisfy the BBN 
bound only if 77™ and 77™ have opposite signs and the total asymmetry is not very different 
from zero. A case with r] u = was the main example presented in [21]. Here we choose instead 
a benchmark case with non-zero negative total asymmetry, 77^ = —0.41, and 77™ = 0.82. The 
evolution of the flavor asymmetries is found from the numerical solution of the EOMs and is 
shown in Figure 1. Other choices of the initial asymmetries will lead to different final values, 
but the overall behavior of the evolution is similar to the case shown here. 

One can see in Figure 1 the effect of flavor oscillations on the evolution of neutrino 
asymmetries with the universe temperature and the dependence on the value of #13. If this 
mixing angle is close to the present experimental upper bound, flavor oscillations are effective 
around T 7 ~ 8 MeV (inner red lines) when neutrinos are still in good thermal contact with 
the ambient plasma. 1 The neutrino distributions evolve keeping an equilibrium form, and 
the total asymmetry is almost equally distributed among the three flavors. However, the 
final value of the electron neutrino asymmetry is too different from zero and this case is not 
allowed by BBN, as we will see in Section 3. Instead, for #13 = flavor oscillations begin 
only at T 7 < 3 MeV, when weak interactions are not frequent enough to keep the neutrino 

1 Here we consider only the case of normal hierarchy, Am^tm > 0. If we choose an inverted hierarchy the 
results are very similar, except that equilibrium among the flavor asymmetries is reached slightly earlier. 
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Figure 1. Evolution of the flavor neutrino asymmetries when 77^ = —0.41, and 77" 1 = 0.82. The solid 
curves correspond to vanishing t9 13 (outer black lines) and sin 2 6*13 = 0.04 (inner red lines). The total 
neutrino asymmetry is constant and equal to three times the value shown (blue dotted line). 



spectra in equilibrium. There is no full equipartition of the total r\ v among the three flavors, 
although collisions lead to final values of the flavor asymmetries closer to r/ u /3 than those 
expected in an MSW matter neutrino conversion. In this case the final value of the electron 
neutrino asymmetry is close enough to zero at the onset of BBN to lie in the favored region. 
If sin 2 #i3 < 10~ 3 , the outcome is very close to the case of vanishing #13. 

In Figure 2 we show the final energy spectra of relic electron neutrinos and antineutrinos 
in arbitrary units for the same case of Figure 1 with vanishing #13. The upper (lower) solid line 
stands for the spectra of electron neutrinos (antineutrinos) calculated numerically, while the 
corresponding dotted lines are described by a Fermi/Dirac distribution just characterized by 
the same effective value of the electron neutrino degeneracy parameter (as used in all analyses 
before [21]). Both cases lead to the same value of the electron neutrino asymmetry but the 
real calculation shows that an excess of radiation in neutrinos remains. 

In Figure 3 we show the evolution of the ratio of neutrino to photon energy densities, 
Pv/p-y, properly normalized so that it corresponds to N e g at early and late times as in [21]. 
The fast drop of Pu/p-y at T ~ 0.2 MeV represents photon heating by e + e~ annihilations. 
The case without asymmetries (dotted line) ends at late times at iV e g- = 3.046 instead of 3 
because of residual neutrino heating [20]. We also show (solid lines) the evolution for our 
main example, where initially N e g = 4.16 for our choice of neutrino asymmetries. One can 
see that as soon as oscillations become effective reducing the flavor asymmetries, the excess 
of entropy is transferred from neutrinos to the electromagnetic plasma, cooling the former 
and heating the latter, but this process is only very effective for large values of #13. While the 
final iV e ff is 3.1 for sin 2 #13 = 0.04, for negligible #13 a significant deviation from equilibrium 
survives and leads to a final enhanced value of iV e g = 3.3. 
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Figure 2. The final energy spectra of relic electron neutrinos and antineutrinos in arbitrary units 
for the same case of Figure 1 with vanishing $13. Upper (lower) solid line stands for electron neutrino 
(antineutrino) calculated numerically (label "real"). Upper (lower) dotted line stands for electron neu- 
trino (antineutrino) described by a Fermi/Dirac distribution just characterized by the same effective 
value of the electron neutrino degeneracy parameter. 

3 Results from BBN: constraints on total and electron neutrino asymme- 
tries 

As well known, BBN depends upon neutrino distribution functions in two ways. First of all, 
electron neutrinos and antineutrinos enter directly in the charge current weak processes which 
rule the neutron/proton chemical equilibrium. A change in the effective temperature of the 
distribution can shift the neutron/proton ratio freeze out temperature and thus modifies the 
primordial 4 He abundance. Similarly, a non-zero v e — v e asymmetry also changes chemical 
equilibrium towards a larger or smaller neutron fraction for negative or positive values of £ Ue , 
respectively. Furthermore, lepton asymmetries in all flavors translate into a positive extra 
contribution to the neutrino energy density, speeding up the expansion rate given by the 
Hubble parameter. 

Differently than in the approximated standard treatments, where both neutrino asym- 
metries and the extra contribution to N e g due to £ Va , see eq. (1.3), are considered as constant 
parameters, in the present analysis we exactly follow the evolution of the neutrino distribu- 
tion versus the photon temperature T~, which is our evolution parameter. 2 To this end we 
have changed the public numerical code PArthENoPE [22, 24] as follows. For any given initial 
values (at T 7 = 10 MeV) for the total neutrino asymmetry r\ v = Y^ a r] Ua , unchanged by fla- 
vor oscillations, and electron neutrino asymmetry rj™ we obtain, as described in the previous 

2 For another example of BBN calculations with arbitrarily-specified, time-dependent neutrino and antineu- 
trino distribution functions, see Ref. [31]. 
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Figure 3. Evolution of the neutrino energy density for the same case as in Figure 1. The vertical 
axis is marked with N c g , left before e + e~ annihilation, right afterwards. The solid curves correspond 
to vanishing #13 (upper black line) and sin 2 813 = 0.04 (lower red line). The case without asymmetries 
is shown for comparison (blue dotted line). 

section, the time dependent neutrino distributions. The latter are then fitted in terms of 
Fermi-Dirac functions with the two evolving parameters T Ua (T^) and ^ Ua (T^). An example 
of their evolution is shown in Figure 4 for the same choice of initial asymmetries as in the 
case described in Section 2. Weak rates are then averaged over the corresponding electron 
(anti)neutrino distribution. The Hubble parameter is also modified to account for the actual 
evolution of total neutrino energy density. 

The final abundances of both the ratio 2 H/H and the 4 He mass fraction, Y p , are nu- 
merically computed as a function of the input parameters r) v and 77™ and compared with the 
corresponding experimental determinations. The baryon density parameter has been set to 
the value determined by the 7-year WMAP result, Q. b h 2 = 0.02260 ±0.00053 (68% C.L.) [5]. 3 

To get confidence intervals for r] v and 7/™, one can construct the likelihood function 

£(^,<)ocexp(-x 2 (^,C)/2) , (3.1) 

with 

X 2 (M£) = £PQ(^,<) -^WyfarfyiXjfarfr) - X° bs ] . (3.2) 

ij 

The proportionality constant can be obtained by requiring normalization to unity, and 
Wij {r\ v , 77™) denotes the inverse covariance matrix [23], 

Wij(r]u, C e ) = Wij + °texp 6 ij + ° 2 ij,other\~ l . ( 3 - 3 ) 
3 The allowed region of fli, in extended cosmological models with free N e g does not differ significantly. 
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Figure 4. Evolution of the effective comoving temperatures and degeneracy parameters of electron 
(solid lines) and muon or tau (dashed lines) neutrinos for the same case as in Figure 1. Both the 
case of vanishing 6> 13 (thick black lines) and sin 2 6* 13 = 0.04 (thin red lines) are shown. The effective 
temperature for the case without asymmetries is shown in the upper panel for comparison (blue dotted 
line). 

where er^- and 

&i,exp represent the nuclear rate uncertainties and experimental uncertainties of 
nuclide abundance X{ , respectively [23] , while by o\- other we denote the propagated squared 
error matrix due to all other input parameter uncertainties (r n , Gn, Qbh 2 , . . . ). In our case 
we consider in eq. (3.2) as Xj the quantities 2 H/H and Y p only. 

Let us now briefly discuss the set of data we have used in our study. The 2 H/H number 
density is obtained by averaging seven determinations obtained in different Quasar Absorp- 
tion Systems, as in [11] 

2 H/H = (2.87 ±0.22) x 1(T 5 , (3.4) 

where the quadratic error has been enlarged by the value of the reduced x 2 to account for 
the dispersion of measurements for this dataset, Xmin/6 = 3.6, see [11]. 

For the 4 He mass fraction we consider two different determinations. One is the result 
of the data collection analysis performed in [11], 

Y p = 0.250 ±0.003. (3.5) 
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Figure 5. The 95% C.L. contours from our BBN analysis in the r\ v — 77™ plane for #13 = (left) 
and sin 2 = 0.04 (right). The two contours correspond to the different choices for the primordial 
4 He abundances of eqs. (3.5) (blue) and (3.8) (purple). The (red) dot-dashed line is the set of values 
of r\ v and 77™ which, due to flavor oscillations, evolve towards a vanishing final value of electron 
neutrino asymmetry 77^". We also report as dashed lines the iso-contours for different values of N e g, 
the effective number of neutrinos after e + e~ annihilation stage. 



More recently, new studies of metal poor H II regions have appeared in the literature [7-9]. 
While these groups both agree on a larger central value with respect to the result of eq. 
(3.5), which incidentally seems to pin down a value for N e g > 3 at 2a, a different estimate of 
possible systematic effects which dominate the total uncertainty budget is quoted in [7] and 
[8], with [8] quoting a larger error, of the order of 4%. 

Y p = 0.2565 ± 0.0010(stat.) ± 0.0050(syst.) [7] , (3.6) 
Y p = 0.2561 ± 0.0108 [8] . (3.7) 

Finally, in a recent paper [9] , Markov Chain Monte Carlo method was exploited to determine 
the 4 He abundance, and the uncertainties derived from observations of metal poor nebulae 
finding 

Y p = 0.2573 ±0.0033. (3.8) 

In the following, we will use the two results of eq.s (3.5) and (3.8). While their uncertainties 
are the same, they differ for the central value, actually the smaller and higher of all results 
reported above, a fact which will produce two different bounds on the electron neutrino 
asymmetry. 

The 95% C.L. contours for the total asymmetry r] u and the initial value of the electron 
neutrino parameter 77™ are shown in Figure 5 for the adopted determinations of 2 H and 4 He 
and for two different choices of #13. In both cases the contours are close to and aligned along 
the red dot-dashed line which represents the set of initial values for the asymmetries which 
eventually evolve toward a vanishing final electron neutrino asymmetry, ff^ ~ 0, which is 
preferred by 4 He data. We recall that 4 He is strongly changed if neutron/proton chemical 



- 10 - 



-1.0 -0.5 0.0 0.5 1.0 

1v 



Figure 6. Bounds in the rf^ - r\ v plane for each nuclear yield. Areas between the lines correspond 
to 95% C.L. regions singled out by the 4 He mass fraction (solid lines) of eq. (3.5), and Deuterium 
(dashed lines) as in eq. (3.4). 
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Figure 7. Same results as in Figure 5 in the plane of initial flavor degeneracy parameters £™ and 

tin 
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equilibrium is shifted by a large value of v e — v e asymmetry around the freezing of weak rates 
(T<y ~ 0.8 MeV). For large #13, oscillations efficiently mix all neutrino flavors and at BBN 
Vv a ~ 77^/3, so the bound on r/ u is quite stringent, —0.1 < T] v < 0.1, if we adopt the value of 
eq. (3.5) for Y p . Instead, for the choice of eq. (3.8) we find —0.20 < rj u ^ 0, i.e. a larger value 
for Y p singles out slightly negative values for rj u (and ry^ n ), since the theoretical prediction 
for Y p grows in this case as the neutron-proton chemical equilibrium shift towards a larger 
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Figure 8. Same results as in Figure 5 in the plane of ?7^ n and the difference rjf, n — rjf, n , where the 
superscript indicates that the asymmetries are evaluated at the onset of BBN, T 7 ~ f MeV, after 
flavor oscillations shuffled the initial ijjf . 





013 = 


sin 2 0i3 = 0.04 


Y p = 0.250 ± 0.003 


-0.66 < rj u < 0.63 


-0.13 <T) V < 0.07 


Y p = 0.2573 ± 0.0033 


-0.71 < r) v < 0.56 


-0.20 < i\u < 0.02 



Table 1. BBN bounds on the initial total neutrino asymmetry at 95% C.L. 



neutron fraction at freeze-out . On the other hand, for a vanishing $13 the contours for r\ v and 
jfy show a clear anticorrelation, and even values of order unity for both parameters are still 
compatible with BBN. The allowed regions of the total neutrino asymmetry are summarized 
in Table 1. 

We stress that for any value of $13, the data on primordial deuterium, eq. (3.4), is crucial 
for closing the allowed region that the 4 He bound fixes along the 77^° ~ line. In fact, though 
2 H is less sensitive than 4 He to neutrino asymmetries and effective temperature which enter 
the Universe expansion rate, see e.g. [11], yet including it in the analysis breaks the degener- 
acy between 77™ and r\ v which is present when only 4 He is used. This can be read from Figure 
6 where the 95% C.L. in the 77™ - r\ v plane are shown for 4 He and 2 H separately, for the case 
013 = 0. The solid lines bound the region of the plane compatible with the 4 He measurement 
as in eq. (3.5), whereas the dashed contours correspond to Deuterium observation, see eq. 
(3.4). The different shape of these two regions is due to the different dependence of nuclide 
abundances on 77J," and r] u , thus their combination breaks the degeneracy and leads to a close 
contour as shown in Figure 5. 

It is also interesting to report our results in terms of other variables, as in Figures 7 
and 8. In the first case, the BBN contours are shown in the plane of initial flavor degeneracy 
parameters £JJ* and £Jf while, in Figure 8, we consider a new pair of variables: the electron 
neutrino asymmetry at the onset of BBN 7y^ n , and the difference 77^" — 77^, which in the 
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standard analysis is usually assumed to be vanishing. One can see from this figure that, 
while the (95% C.L.) bound on ry™ 1 is independent of the value of #13 

-0.03 < rf^ < 0.02 , for Y p = 0.250 ± 0.003 , (3.9) 
-0.05 < rf^ < 0.003 , for Y p = 0.2573 ± 0.0033 , (3.10) 

the difference between the final v e and v x asymmetries strongly depends upon this yet un- 
known mixing angle, as expected. In fact, for large #13 we recover the standard result, 
r]~ a ~ r]„, due to efficient mixing by oscillations and collisions, while for #13 = the two 
asymmetries can be different. The following ranges are in good agreement with BBN data 
at 95 % C.L., independently of the adopted value for Y p 

- 0.3 < rftl - rt e < 0.3 , sin 2 9 13 = , (3.11) 
-0.04 < nil - rfll < 0.04 , sin 2 13 = 0.04 . (3.12) 

We conclude that, in particular for #13 = oscillations lead to quite different neutrino 
asymmetries in e and fi/r flavors, still being in good agreement with BBN, differently than 
what was previously assumed in the literature, see [1-4], and [11-19] for BBN analyses. 

In Figures 5-8 we also plot iso-contours for the value of the effective number of neutrinos, 
iV e ff , evaluated after e + e~ annihilations. For large #13 BBN data bound N e g to be very close 
to the standard value 3.046, since all asymmetries should be very small in this case and flavor 
oscillations modify the neutrino distributions while neutrinos are still strongly coupled to the 
electromagnetic bath. Therefore, we do not expect non-thermal features in the neutrino 
spectra in this case, since scatterings and pair processes allow for an efficient transfer of any 
entropy excess. On the other hand, for vanishing #13, larger values of N e s are still compatible 
with BBN data, up to values of the order of 3.4 at 95% C.L. The dependence of the largest 
achievable value of N e g on the value of the mixing angle #13, obtained by spanning in the 
asymmetry parameter plane the region compatible with BBN, is reported in Figure 9. 

It is worth noticing that the final values of N e g, in particular for large final asymme- 
tries of is x , are also slightly larger than the N e s that one would obtain using the equilibrium 
expression of eq. (1.2). For example, if we take 7/JJ n = and 77JJ 11 = 0.3, a point on the bound- 
ary of the BBN contours (see Figure 8) and compute the corresponding effective chemical 
potentials, using eq. (1.3) one gets N c s = 3.2, while the actual value is larger, N e s = 3.4, 
a signal that in this case the interplay of solar-like oscillations and neutrino freeze-out has 
produced indeed, a mild non-thermal distortion in neutrino distributions. 

Finally, we notice that our BBN results correspond to a minimal scenario with primordial 
neutrino asymmetries, as we do not consider a possible extra contribution to iV e g coming from 
relativistic degrees of freedom other than standard active neutrinos. Their effect is known 
to produce looser bounds on neutrino asymmetries, as they speed up expansion and thus, 
can compensate the effect of a positive v e — v e asymmetry. We have explicitly checked that, 
for some choices of primordial asymmetries, the addition of extra radiation does not modify 
the evolution of flavor neutrino asymmetries. Of course, in such a case the contribution to 
the energy density of the additional relativistic degrees of freedom adds up to the surviving 
excess to N e g arising from neutrino asymmetries. 

4 Conclusions 

In this paper we have calculated the evolution of neutrinos in the early universe with initial 
flavor asymmetries, taking into account the combined effect of collisions and oscillations. Our 
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Figure 9. Largest values of iV e ff from primordial neutrino asymmetries compatible with BBN, as a 
function of #13. 



numerical results for the neutrino momentum distributions, which can develop non-thermal 
features, were used to find the primordial production of light elements, employing a modified 
version of the PArthENoPE BBN code. Comparing with the recent data on primordial 2 H and 
4 He abundances, we have found the allowed ranges for both the total asymmetry and the 
initial asymmetry in electronic flavor. These BBN limits mostly depend on the value of $13, 
the only unknown mixing angle of neutrinos. As can be seen from Figure 9, for sin 2 #13 < 10~ 3 
this implies an effective number of neutrinos bound to N^g < 3.4, independently of which 
of the experimental 4 He mass fraction of eq.s (3.5) and (3.8) is adopted, whereas for larger 
values of #13 the effective number of neutrinos is closely bound to the standard value 3.046 
obtained for vanishing asymmetries. 

In the near future it will be possible to improve our BBN constraints on the lepton 
number of the universe. On one hand, thanks to the better sensitivity of new neutrino 
experiments, either long-baseline or reactor, we expect to have a very stringent bound on #13 
or eventually its measurement (see e.g. [32]). Such results would lead to a more restrictive 
BBN analysis on primordial asymmetries. On the other hand, in the next couple of years 
data on the anisotropies of the cosmic microwave background from the Planck satellite [33] 
will largely reduce the allowed range of iV e g, since the forecast sensitivity is of the order of 
0.4 at 2a [34-36]. Indeed, suppose that Planck data confirm a value of iV e g > 3. An excess of 
radiation up to 0.4 — 0.5 could imply the presence of a significative degeneracy in the neutrino 
sector, but this would depend on the result of #13 measurement. A vanishing #13 would allow 
such possibility, but a measured #13 in the next generation of experiments would imply the 
presence of extra degrees of freedom others than active neutrinos. In case of a much larger 
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result by Planck, namely N e g > 4, it would be impossible to explain such a result in terms of 
primordial neutrino asymmetries only, and alternative cosmological scenarios with additional 
relativistic species, such as sterile neutrinos (see for instance [37]) would be strongly favored. 
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